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Abstract. Analyses of recent helioseismic data indicate 
that the dynamical regimes at the base of the convec- 
tion zone can be different from those observed at the top, 
having either significantly shorter periods or non-periodic 
behaviour. 

Recently spatiotemporal fragmentation/bifurcation 
has been proposed as a dynamical mechanism to account 
for the multi-mode behaviour that is possibly observed 
in the solar convection zone, without requiring separate 
physical mechanisms with different time scales at differ- 
ent depths. 

Here we study the robustness of this mechanism with 
respect to changes to the zero order rotation profile, moti- 
vated by the uncertainties of and differences between the 
various reductions of the helioseimic data. We find that 
spatiotemporal fragmentation is a common feature of the 
reductions investigated. 
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1. Introduction 

Recent analyses of the helioseismic data, from both the 
Michelson Doppler Imager (MDI) instrument on board the 
SOHO spacecraft (Toomre et al. 2000) and the Global Os- 
cillation Network Group (GONG) project (Antia & Basu 
2000), have provided strong evidence that the previously 
observed torsional oscillations (e.g. Howard & LaBonte 
1980; Snodgrass, Howard & Webster 1985; Kosovichev & 
Schou 1997; Schou et al. 1998), with periods of about 11 
years, penetrate into the convection zone to depths of at 
least 10 percent in radius. 

Further studies of these data have produced interest- 
ing, but rather inconsistent results. In particular, Howe et 
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al. (2000b) find evidence for the presence of similar oscil- 
lations near the tachocline situated close to the bottom of 
the convection zone, but with a markedly shorter period 
of about 1.3 years, whereas Antia & Basu (2000) do not 
find such oscillations (whilst allowing the possibility that 
more irregular variations may be present). Whatever the 
true dynamical behaviour at these lower levels turns out 
to be, both these results indicate that the variations in 
the differential rotation can have different dynamical be- 
haviour at the top and the bottom of the solar convection 
zone, with the oscillations at the bottom having either 
significantly shorter periods or non-periodic behaviour. 

Recently, spatiotemporal fragmentation/bifurcation 
has been proposed as a dynamical mechanism to account 
for the possible multi-mode behaviour in different parts 
of the solar convection zone (Covas et al. 2000a, hereafter 
CTMl). Evidence for this mechanism was found in the 
context of a two dimensional axisymmetric mean field 
dynamo model in a spherical shell, with a semi-open 
outer boundary condition, in which the only nonlinearity 
is the action of the azimuthal component of the Lorentz 
force of the dynamo generated magnetic field on the 
solar angular velocity. The underlying zero order angular 
velocity was chosen to be consistent with the recent 
helioseismic data. 

Our initial study (CTMl) suffered from two major 
shortcomings. The dynamo model employed inevitably 
contained major approximations, and the zeroth order ro- 
tation profile used is bound to include uncertainties, given 
the uncertainties in the inversion techniques, as well as the 
short extent in time over which data sets are so far avail- 
able. 

To address the first shortcoming, an extensive study 
was made of the persistence of the fragmentation following 
plausible changes in the details of the dynamo model. This 
study showed that the mechanism remains surprisingly 
robust in presence of a variety of rather severe changes to 
the model (Covas et al. 2001, hereafter CTM2). 

Regarding the second shortcoming, CTM2 also in- 
cluded a preliminary study, allowing the zero order rota- 
tional profile to be given by either the inversion of the MDI 
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data reported by Toomre et al (2000) or of the correspond- 
ing GONG data (Howe et al 2000b). This also showed the 
mechanism to be robust with respect to such changes to 
the rotation profile. 

Given the significant uncertainties that still remain in 
the helioseismic measurements, especially given the lim- 
ited length of the data so far available, it is crucial to 
study the robustness of the spatiotemporal fragmentation 
mechanism to a wider range of the rotation profiles con- 
sistent with the observations. 

2. Zero order rotation profile 

New measurements of the solar internal rotation with 
SOHO MDI data have been reported recently by 
Vorontsov et al (2001), using novel inversion techniques 
that are quite different from those employed earlier. Re- 
sults of a 2D inversion confirmed the previous findings that 
the torsional oscillations (migrating zonal flows) penetrate 
deep into the convection zone. However, given the limited 
accuracy of the helioseismic data currently available, they 
did not allow the resolution of these oscillations all the way 
down to the bottom of the convection zone. The so-called 
1.5D inversion procedure has also been used, which con- 
sists in measuring separately the consecutive terms fli (r) 
in the expansion 
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where P2i+\ is the (2i -I- l)th Legendre polynomial, r is 
the distance from the centre and d is the co-latitude. 

The 1.5D inversion of Vorontsov et al (2001) was prin- 
cipally targeted at improvement of the spatial resolution 
of the rotation profile near and just below the base of 
the convection zone, in the tachocline region — the region 
where the dynamo mechanism is most sensitive to the ro- 
tation profile, since the largest gradients of the angular 
velocity appear to be localized there. 

An important feature of the solar internal rotation as 
seen in the helioseismic measurements is that the rapid 
spatial variation of the angular velocity in the tachocline 
is governed principally by the first two or three terms 
in the expansion (1); the contribution of higher terms 
appears to be an order of magnitude smaller, being es- 
sentially undetectable in the helioseismic data available 
at present. Another important feature is that this domi- 
nant part of the variation of the angular velocity in the 
tachocline with depth and latitude, captured by the first 
few terms in the expansion (1), does not show any changes 
with solar activity (Vorontsov et al 2001). Taken together, 
these two features suggest considering the torsional oscil- 
lations as a small-amplitude, rapidly- varying (in Q\ time- 
dependent component of the solar rotation, imposed on a 
large-amplitude, slowly- varying (in d\ nearly-stationary 
component described by the first terms of the expansion 
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These considerations show that at least a dominant 
part of the time-independent (background) component of 
solar rotation can be evaluated by a simple truncation of 
the expansion (1) to the first two or three terms, mak- 
ing the truncated expansion a natural candidate for the 
zero-order rotation profile in the dynamo modelling. This 
choice might be more consistent than using a time average 
of the inverted rotation profiles (as in CTM2), since the 
observational data used in the inversions do not yet cover 
enough (ca. 5 yr) of a complete solar cycle to allow proper 
averaging. With the zero-order profile derived from the 
time-independent part of the data we can also eliminate 
any possibility that torsional oscillations produced by the 
dynamo model originate from the time varying zonal flows 
already present (in time-averaged form) in the background 
rotation profile. 

Thus we here describe a comparative study of the ef- 
fects of employing zero order rotation profiles, obtained 
by taking the first two and the first three terms of the ex- 
pansion (|l|). These profiles, inferred from the 360d SOHO 
MDI data set, are shown in Fig. [|. 

3. The model 

In order to test the sensitivity of the spatiotemporal frag- 
mentation with respect to the uncertainties in the rotation 
profile, we used the following dynamo model, which was 
also used in CTM1,2. We assume that the gross features 
of the large scale solar magnetic field can be described by 
a mean field dynamo model, with the standard equation 



SB 

dt 



= V X (u X B + qB - r;V X B). 



(2) 



Here u = w0 — ^V?], the term proportional to V77 repre- 
sents the effects of turbulent diamagnctism, and the ve- 
locity field is taken to be of the form w = + v', where 
vq — r^o'^sinfl, r^o is a prescribed underlying rotation law 
and the component v' satisfies 



dv' _ (V X B) X B 
dt \iQp 



(3) 



where is the operator 
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and hq is the induction constant. The sole nonlin- 
earity in the dynamo equation is the feedback of the az- 
imuthal component of the Lorentz force (Eq. (^), which 
modifies only slightly the underlying imposed rotation 
law, but nevertheless limits the magnetic fields at finite 
amplitude. The assumption of axisymmetry allows the 
field B to be split simply into toroidal and poloidal parts, 
B = Bt Bp = V X and Eq. (|) then yields 

two scalar equations for A and B. Nondimensionalizing in 
terms of the solar radius R and time R^/rjQ, where rjo is 
the maximum value of i], and putting O = a = aoa, 
V = VoV, B = -BqB and v' = n*Rv', results in a sys- 
tem of equations for A, B and v' . The dynamo parameters 
are the two magnetic Reynolds numbers Ra — a^R/riQ 
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Fig. 1. Isolines of the angular velocity of the solar rota- 
tion, obtained from the 1.5-D inversion of the 360d SOHO 
MDI data set by Vorontsov et al (2001), given by the first 
two terms (top panel) and the first three terms (bottom 
panel) of the expansion (|l|) . Contours are labelled in units 
of nHz. 

and = n*R'^/rio, and the turbulent Prandtl number 
Pr — i^a/vo- ^* is the solar surface equatorial angular ve- 
locity and 77 — rj/rjQ. Thus vq and 770 are the turbulent 
magnetic diffusivity and viscosity respectively, Ruj is fixed 
when 770 is determined (see Sect. but the value of Ra is 
more uncertain. The density p is assumed to be uniform. 

Given the ill-determined nature of boundary condi- 
tions in astrophysical settings, we make physically moti- 
vated choices following CTM1,2. At inner boundary we 
chose B = 0, ensuring angular momentum conservation, 
and an overshoot-type condition on Bp (cf. Moss & 
Brooke 2000). At the outer boundary, we used an open 
boundary condition dB/dr = on i? and used vacuum 
boundary conditions for Bp (see Kitchatinov et al. (2000) 
and CTM1,2 for motivation). 

Equations (|^) and (|3|) were solved using the code de- 
scribed in Moss & Brooke (2000) (where more details 
are given; see also Covas et al. 2000a) together with the 
above boundary conditions, over the range < r < 1, 
< < TT. We set ro = 0.64; with the solar convection 



zone proper being thought to occupy the region r > 0.7, 
the region < ^ 0.7 can be thought of as an overshoot 
region/tachocline. We note that there is some evidence 
that the tachocline may be rather narrower than can be 
directly resolved by the usual inversion techniques. The 
uncertainties associated with this possibility are difficult 
to assess, but we do note that the appearance of fragmen- 
tation appears to be quite robust with respect to changes 
in the model (CTM2). In the following simulations we used 
a mesh resolution of 61 x 101 points, uniformly distributed 
in radius and latitude respectively. 

In this investigation, we took the zero order rotation 
profile fi(r, 9) in the region 0.64 < r < 1 to be given 
successively by the first two and the first three terms of 
the series (|]) given above. 

For the alpha-effect, we took a — ar{r)f{d), where 
f{9) = sin"* cos (cf. Riidiger & Brandenburg 1995) and 
Ur — 1 for 0.7 < r < 0.8 with cubic interpolation to zero 
at r = ro and r — 1, with the convention that a,. > and 
Ra < 0. Also, in order to take into account the likely de- 
crease in the turbulent diffusion coefficient rj in the over- 
shoot region, we allowed a simple linear decrease from 
77 = 1 at r = 0.8 to 77 = 0.5 in r < 0.7. We note, however, 
that as was shown in CTM2, these latter details of the 
model are unlikely to remove the possibility of spatiotem- 
poral fragmentations. 



4. Results 

We calibrated our model so that near marginal excita- 
tion the cycle period was about 22 years. This determined 
R^ = 44000, corresponding to 770 « 3.4 x 10^^ cm^ sec~^, 
given the known values of fl* and R. With the rotation 
law given by the first three terms of (1), the first solutions 
to be excited in the linear theory are limit cycles with 
odd (dipolar) parity with respect to the equator, with 
marginal dynamo number R^ ~ —3.76. The even par- 
ity (quadrupolar) solutions are also excited at a similar 
marginal dynamo number of Ra ~ —3.95. These values 
did not change significantly when the first two terms of 
expansion (1) were used. It is arguable that the turbu- 
lent Prandtl number is of the order of unity, and we set 
Pr — 1. For the parameter range that we investigated, the 
even parity solutions can be nonlinearly stable. Given that 
the Sun is observed to be close to an odd (dipolar) par- 
ity state, and that previous experience shows that small 
changes in the physical model can cause a change between 
odd and even parities in the stable nonlinear solution, we 
chose to impose dipolar parity on our solutions. 

With these parameter values, each of the rotation pro- 
files considered here produced butterfly diagrams which 
are in qualitative agreement with the observations. The 
model also successfully produced torsional oscillations 
that penetrate into the convection zone, in all cases con- 
sidered, similar to those deduced from recent helioseismic 
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data (Howe et al. 2000a, Vorontsov et al. 2001), and stud- 
ied in Covas et al. (2000b) 



plotted in Fig. g the radius of the top of the fragmentation 
region. This qualitative similarity of the fragmentation in 
all the cases is further evidence of robustness. 
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Fig. 2. Radial contours of the angular velocity residuals 
Sn as a function of time for a cut at latitude 30°, with 
the zero order rotation profile obtained by taking the first 
two (top panel) and the first three terms (bottom panel) 
in (|^). Parameter values are Ra = —20.0, P,. = 1.0 and 
= 44000. 

Both of the rotation profiles considered here produced 
spatiotemporal fragmentation. As examples, the top and 
the bottom panels of Fig. ||show the radial contours of the 
angular velocity residuals Sfl as a function of time for a cut 
at latitude 30°, with the zero order rotation profile given 
by the first two and the first three terms in the expression 
(|l|), respectively. 

To test further this robustness, we also used func- 
tions fl2i+i{r)^i = 0,1,2, obtained from an inversion 
with slightly different regularisation parameters and in 
a somewhat different frequency range of the input data 
(Vorontsov et al 2000). We found no qualitative changes 
in the spatiotemporal fragmentation. 

Finally, as a comparison with our previous results ob- 
tained by employing the GONG and MDI data, we have 
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Fig. 3. The height of the top of the fragmentation region 
for different zero order rotation profiles flo and different 
dynamo parameters Ra- 



5. Discussion 

We have shown that spatiotemporal fragmentation per- 
sists even when only the lowest order terms in the Legen- 
dre expansion of the inversion for the solar rotation law 
are used. This is important, as these terms appear to be 
approximately time independent (over the scale of a so- 
lar cycle) , and also dominate the expansion. Indeed, these 
terms really represent the rotation law that should be in- 
put to the dynamo code - in principle solution of Eq. ^ 
should then determine the variations of fl over a cycle. 
Thus, these results allay worries that the inconsistency of 
averaging the observed rotation law over only a part of a 
solar cycle might be giving misleading results. 

The results presented here, together with those of 
CTMl and CTM2, show spatiotemporal fragmentation to 
be a very persistent feature of our simulations. This find- 
ing suggests that nonlinear dynamo modelling has a real 
potential predictive power, which could lead to a better 
understanding of the solar dynamo when longer data sets 
as well as improved analyses of them become available. 
This in turn, in conjunction with the expected improve- 
ments in the helioseismological data, could lead to an im- 
proved determination of the properties of the solar inte- 
rior, including the nature of the tachocline. 
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